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ABSTRACT 


False Sonar targets present a serious unpredicted problem to U.S. 
Navy ASW units. It is гс на that planning and operations could be 
enhanced by a forecasting capability for whale distribution. As a 
possible solution to this problem, a modified form of the "Transect 
Method of population estimation" is applied to whaling data to calculate 
probable numbers of false targets per 1000 nautical miles of steaming 
with a 1000 yard sonar range. Japanese and Russian whale fishery data 
are analyzed by the "q" and Expected Catch methods of population dynamics 
to obtain two independent estimates of the populations of fin, sei and 
sperm whales. The mean of the two estimates is applied to the equation 
along with a term for assumed ideal sonar conditions. The data is 
calculated by ten degree square of latitude and longitude, north of 30°N, 
and presented on Fleet Numerical Weather Central polar stereographic 
charts for the months April through December. The number of false targets 
attributable to fin, sei and sperm whales alone range from 1 to 63 south 


of the Aleutian Islands and 1 to 30 off Honshu, Japan. 
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1. INTRODUCTION 


To the ASW unit, one of the most serious problems encountered is the 
vast amount of time wasted prosecuting spurious or false sonar contacts. 
Detection of a contact requires the detachment of a ship or ships from 
the task group to investigate the threat; thereby diminishing the overall 
capability of the screen in protecting the main body. Of the many 
sources of false sonar contacts, the most troublesome and most difficult 
to distinguish from actual submarine contacts are whales. It is, there- 
fore, imperative that methods be developed to reliably predict areas of 
high whale density which will serve to alert sonar operators. 

Since the resumption of the Japanese whale fishery following World 
War II, a great deal ee has been generated in determining the 
areas of whale abundance and seasonal fluctuations in whaling grounds, 
resulting in a research effort in the sub-arctic North Pacific Ocean. 

To date, research has proceeded in three major areas of interest: 
geographical distribution of the major whale species, investigation of 
biological and behavioral characteristics, and examination of the 
oceanographic conditions of the whaling grounds. Each of the three 

is essential in attempting to determine the factors which control whale 
migration and defining favorable whaling regions. 

The first major work on whale migrations (Kellog 1929) described the 
primary paths of migration for blue (Balaenoptera musculus), sei 
(B. borealis) and fin (B. physalus) whales in the Pacific and Atlantic 
Oceans, in general terms, based on all sighting and catch information 


available at that time. Whaling grounds developed over the years and 





were illustrated by Townsend (1935) by plotting catches according to 
logbook records of American whaling ships for a period of one hundred 
fifty nine years (1761-1920). Other researchers have concentrated on 
more limited areas of the North Pacific in response to the whaling 
interests of their own country's territorial waters (Omura and Nemoto 
mesos Clarke 1957; Sleptsov 1961; Tarasevich 1965, 1970; Nishiwaki and 
Kasuya 1970; Kasuya 1971). General distribution charts of cetaceans in 
the North Pacific were developed by Nishiwaki (1967). Experimental 
investigations were also conducted (and continue) by means of whale 
marking (Shevill and Watkins 1966; Ivashin and Rovnin 1966). 

Feeding habits and food preferences were investigated by observation 
and investigation of stomach contents (Nishimoto, Tozawa and Kawakami 
1952; Nemoto 1957; Nemoto 1959; Nemoto 1962; Nemoto and Kasu 1965) which 
were correlated with migration cycles of baleen whales (Nemoto 1959). 
Schooling behavior in the feeding areas was examined (Nemoto 1964) and 
compared with patterns of sexual segregation (Ohsumi 1966). 

Consideration of biological parameters has mainly been confined to 
development of reliable methods of age determination by ear-plug lamina- 
tion count (Nishiwaki 1952; Nishiwaki and Masaharu 1952; Nishiwaki et al 
1957; Nishiwaki 1957; Nishiwaki 1958; Ichihara 1969, Ichihara 1963, 
Ichihara 1964; Ichihara 1966) and dentine layer count in the case of 
sperm whales (Nishiwaki 1958; Ohsumi 1963; Gambell and Grzegorzewska 
1967). Complimentary studies have focused attention on the determination 
of average body proportions and growth rates (Mathews 1938; Omura 1950; 
Nishiwaki and Ohe 1951; Nishiwaki and Hibiya 1951; Nishiwaki et al 1952; 
Ohno 1952; Omura and Fujino 1954; Nishiwaki 1954; Fujino 1954; Omura and 


Nemoto 1955; Fujino 1955; Nishiwaki 1956; Ohsumi 1958; Nasu and Masaki 1970). 





Uda (1954), Uda and Dairokuno (1957), Uda and Suzuki (1958), Nasu 
(1960), Uda (1962), Nasu (1963), Nasu (1966) have conducted oceanographic 
dies with the objective of determining the exact physical parameters 
which define favorable whaling grounds and correlating these with meter- 
ological conditions (Uda 1956) e.g. fluctuations in sea surface temperature, 
salinity and nutrients. 

.Leapley and Levenson (1969) have published the only study to date 
eveltvating the false target threat. This ee contours of expected 
numbers of false targets due to whales in the eastern North Atlantic ocean, 
based on sighting records obtained over a five year period. 

The objective of this thesis is to take advantage of these studies to 
provide information on the number of whales which an ASW ship may expect 
to encounter as false sonar contacts during operations in the North 
Pacific Ocean. Instrumental to this objective is the application of a 
modified form of the "transect method of population estimation" to whale 
fishery data reported by Japanese and Russian whaling companies north of 
30°N. This information may be applied by fleet units operating in the 


North Pacific and by commanders planning operations in these areas. 
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11. METHODS 


A recognized and frequently applied method of population estimation 
for game, birds, and other animals, based on sighting observations, is 
the Transect Method. This method uses the number of animals sighted in 
a known distance, combined with the distance from the track of the 
Sightings to plot contours of relative abundance. 

Mackintosh and Brown (1956) developed a modified version of this 
method for use in determining the fin whale population of the Antarctic 
whaling grounds. It is based on the number of sightings made while 
steaming through a representative distance of an ocean area when no 
estimates are made of the distance of the whales from the track of the 
Sighting vessel. This и may be expressed: 


| A(nni?) n (whales) 
N (whales) - ЖЕТЕ | яв он 


where: 

N = calculated total population of area A 

A = area of ocean considered in square nautical miles, nni 

L = distance steamed through the area in nautical miles, nmi 

D - lookout visibility in nautical miles, nmi 

n - number of whales sighted in distance L 

p = an empirically determined probability factor based upon 

weather conditions, sea state and visibility. 

Consideration of the behavioral characteristics of whales in feeding 

and migration suggested the applicability of this method to evaluation 


> 


of false target threats. Nasu (1966) reports that feeding whales are 
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commonly found at depths less than fifty meters. The exception to this 
would be sperm whales, which occasionally dive to great depths to feed on 
giant squid (Heezen 1975; Okutani and Nemoto 1964). However, in the 
North Pacific, they most frequently feed on schooling fish and therefore 
can also be considered to be mainly distributed in the upper fifty meters. 
Migration also takes place on or near the surface. Thus, although a sonar 
ensonifies a volume of water, in the case of false targets attributable 
to whales we may consider this ensonified region as an area. It is then 
readily apparent that the searching sonar beam is analogous to the eye 
of the lookout on a surface ship and this method may be applied to the 
evaluation of false targets. In this concept, the sonar range is analogous 
to the lookout's visibility. 

Assuming ideal sonar conditions: 


n (number of false targets) - P Lm) x Denm) x N (whales) 
A (nm?) 
where: 
n = expected number of false targets 
L = distance steamed in nautical miles 
D = two times the sonar range in nautical miles 


A = area of the oceanic region 


N = calculated population of area A 


P = probability factor determined by sonar conditions. 
These quantities are illustrated in Fig. 1. 

For the purposes of this study, ideal sonar conditions were assumed 
and, therefore, P was taken to have a value of one; indicating that if 


a whale is present it will be detected. Under actual operating condi- 


tions this criterium is rarely met since the value will seldom be equal 
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Figure 1. Illustration of Terms Used in Modified Transect 
Equation: A = ocean area, D = sonar range, 
L = distance steamed through area A. 





to one, but will be determined by the factors affecting the sonar range 


in the particular area of operation. 


In order to make clearer the techniques used and the sequence of 
operations, the analysis scheme followed in this thesis is presented 


diagramatically in the following flow diagram. 


FLOW DIAGRAM 


RAW DATA 









"q" METHOD [EXPECTED CATCH METHOD] 
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PROGRAM AGEPROP 


PROGRAM WHLPOP 


POPULATION ESTIMATES 





POPULATION ESTIMATES 


MEAN POPULATION ESTEATE 


MODIFIED TRANSECT METHOD 











FALSE TARGET PROGRAM] 


EXPECTED FALSE TARGET DISTRIBUTION 
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The raw data is taken from Russian and Japanese North Pacific fin 
M@alaenoptera physalus), seil (B. borealis) and sperm (Physeter catadon) 
Шише fishery catch results for the years 1966 to 1970. The total catch 
was listed by ten degree square, effort (in catcher days) and length 
frequency distribution. Two independent methods of population dynamics, 
“Іле "а"! method and the expected catch method, ECM, were utilized to 
estimate populations in selected geographical subdivisions shown in Fig. 
Ze For the purposes of this study, Brydes PM were included with sei 
whales since they are distinct from sei whales only by the presence of 
ventral grooves and otherwise have the same physical characteristics. 

The age composition of the catch was determined and the data analyzed 
by the "q' method (Allen 1966), utilizing a computer program written to 
apply this method. This program, designated WHLPOP (Appendix B), yielded 
Hie initial population estimate. 

The second independent estimate of population size was calculated by 
means of the "Expected Catch" method (Allen 1966) using values for natural 
mortality and recruitment! obtained from the reports of the International 
Whaling Commission (IWC 1966-1970). This method was applied using a 
program provided by Dr. K. R. Allen, Fisheries Research Board of Canada, 
Biology Station, Nanaimo, British Columbia. The program was modified to 
conform to the IBM 360 located at the W.R. Church Computer Center, Naval 
Postgraduate School, Monterey, California. 

In accordance with procedures followed by the International Whaling 


Commission (IWC 1968) the mean of the two population estimates was taken 


1 : | 
The addition of whales to the hunted portion of the population by 
growth from younger ages. 
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and the average expected number of whales in each ten degree square, for 
the period 1966-1970, was calculated by means of proportions as expressed 


in equation (1). 


э (1) 


СА 1 


Number of whales in square (М) = 





where: 
М, = estimated mean whale population in the selected area by "q" 


and ECM for the five year period 1966 to 1970. 


C. = average total catch in the ten degree square over the five 
years. 
Ca = average total catch in the selected area 1966-1970. 


Similarly, the average number of whales present per month in each 


ten degree square was calculated by equation (2): 





Number of whales per month (N ) = E X No (2) 
S 
where: 
Е = average total catch in the month 
Co = average total catch im the ten degree square 
N. = population of the square calculated in equation (1). 


Historical records of Japanese whale fishing (Figures 3-5) for the 
years 1945 to 1962 (Nishiwaki, 1966, 1967) were combined with studies of 
migratory cycles and distribution by Nishiwaki and Kasuya (1970) and 
historical catches (Townsend 1955). This was done to determine areas of 
distribution for fin, sei and sperm whales. From these, a composite area 
chart was developed, Figure 6, showing where whales could be encountered 
as false targets. In areas where no historical geographical data was 


= 


available, i.e. Russian fisheries in squares M28-30 and N27-30 (Figure 2), 
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a uniform distribution of the calculated average whale population per 
month was assumed over the entire ten degree square. The number of 
whales per one thousand miles steaming with a 1000 yard sonar range was 
then calculated using the modified form of the Transect Method of popula- 
tion estimation described earlier. This was then plotted on Fleet 


Numerical Weather Central polar sterographic charts (Figures 24-32). 


СИЕ ТИЕ "q" METHOD 

Whaling in the North Pacific produces data which is highly variable, 
Particularly in the measure of catch per unit 2. This variation is 
produced primarily as a result of the smaller effort employed, as compared 
Ine Antarctic fishery, and the exploratory nature of the fishing effort 
to develop new whaling grounds. Methods utilized in population estimation 
must therefore be designed to normalize the data in order to give reliable 
e amates even in the presence of this high variability. One such method 
tithe "q" method. 

The basic assumption of the "а" method is that the ratio between the 
instantaneous fishing mortalitv rate, F, and the total fishing effort, 


X, is constant, i.e. 


q = instantaneous fishing mortality _ CONSENT (3) 
ЕПИР Errort 


Using the pooled data for as many years as possible, an estimate of "q" 
1s obtained and the population computed by dividing the catch per unit 

effort, C/E, by q. In this method a base year is selected in which it 

is assumed that the population structure is determined by natural 


mortality and that recruitment has been constant from year to year. 


2 
The total catch of whales, in numbers, divided by the effort 
required to produce the catch (equation 4). 
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, 


The essential features of this method are: 

CDA Single year group? is followed throughout to avoid the effects 
e recruitment. 

(2) The relative abundance of the selected year or age group is 
compared in each succeeding year with the base year4, where it is assumed 
that the age structure is unaffected by exploitation. Thus, the population 
structure by ages is determined by natural effects with no deaths due to 
fishing. 

(3) Estimates of the total fishing mortality”, for all age groups, 
and the instantaneous fishing mortality of the specified age group are 
Seamed for increasing periods of years. Thus the effects of data 
fluctuations due to seasonal variation are minimized. These estimates 
rc inen divided by the total effort over the same periods to obtain a 
Meameecstimate of "a," table II. 

(4) The effort expended in the whale fishery by different size 
catchers with varying endurance is standardized by multiplying catcher 


days by vessel tonnage and dividing by one thousand, i.e. 


Бро tener days mort catcher tonnage (4) 


The portion of the population consisting of all whales of a 
specified age and older. 


4 Ап arbitrarily selected year in which the population is assumed 
to be affected by natural effects only with no deaths attributable to 
the effects of man. This year is used as a reference. 


The ratio of the number of deaths to the total number of whales in 
the population, summed over the year. 
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The basic equation is: 


Q 
Е = Ш l,r x Qo ,r+t-1 x 01 (5) 
t Qt,r+1-1 Q0,r J 
Wiere: 
F. = total instantaneous fishing mortality from year 1 to year t. 


ET = proportion of the catch in year 1 of whales of age r 
(the first fully recruited=age-group)-and older. 

Q. r+t-1 e proportion ofthe catch in year t of this same year 
eroup (i.e 11 “Іш the base year six year old whales are 
examined then two years later the proportion of whales 
eight years and older must be considered. 

Qo,r+t-1 = proportion from the catch of the previously unexploited 
stock of the same age groups. 


Qo > = proportion in the catch from the previously unexploited 
2 


stock of whales of age r and older. 


Q, = catch per unit effort in year 1. 
9c 


lese quantities are illustrated in table І. 


eatcehspersunitzeitfortzın year t. 


Since the estimate of F obtained in this method refers to the 
decrease from the average population in one season to the average level 


6 


in a later year, the effort X^ may be written by means of integrating by 


the Trapazoid Rule: 


EI 
Х, х пе + 2X. 4 + «| (6) 


The total ene rrort Е, expended from mid-season in one year 
to mid-season in the next. 
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TABLE I 
ILLUSTRATION*OF QUANTITIES USED IN ''q'" METHOD 


Cumulative Proportions in Catch of Whales 
Age (r) and Older (Russian Area V) 


Age 
Group Base Year Leap meat 2 Year 5 Year 4 
66 5 1967 5 1968 5 1969 5 ]970 5 
0 1.0000 1.0000 1.0000 - 1.0000 1.0000 
1 . 9984 1.0000 1.0000 1.0000 -9960 
2 19892 29991 . 9967 . 9969 . 0799 
3 . 8768 200905 .8876 . 6220 . 8556 
4 5227 . 8926 20509 .7604 7566 
5 . 7438 . 8043 - 7256 20557 .6468 
6 .6208Q. .r .68100 257112 . 5260 25225 
7 ‚46720 5024 275 39850 3591 3527 

O. rrt=1 з : t,r+t-1 ; : 

8 SÓ . 4125 . 3003 22121 12769 
9 · 2880 22250 22021 · 1888 · 1862 
10 „2105 „2/54 ‚1414 . 1257 059 
11 . 1357 . 2242 .0812 .0686 . 0861 


Moral] set of data are in appendix. 


*Qutput from Computer Program AGEPROP in Proportion of Each Age 
Group in the Total Catch Converted to Cumulative Distribution 
mem Purposes of Illustration. 
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The error due to random effects on the fishery may be minimized by 
extending the estimate of "q'' from the base year to successive years and 
taking the mean as the value for the area, Table II. 


Thus: 


SE 


q = 7 
=X, (7) 


мар ма 


where T is the total number of years inthe series. 


An estimate of the mid-season recruited population is then given 


by: 


The "q" method may be explained more simply by the block diagram in 
Figure 7. 
FIGURE 7 


Block Diagram Representation of the "а" Method 


BASE YEAR 





M 


natural mortality 


F 


fishing mortality 
The large rectangles represent the total catches in years 1 and 2. 
Considering first the catch in year 2, the small rectangle, 7*, represents 


= 


the proportion of the catch in year 2 which is seven years old and older. 
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It can easily be seen that if there were no fishing or natural mortality 
then the level for the seven year age group would include rectangles F 

and M. Going back one year, the whales which were seven years old and 
older in year two are represented by the six year old and older whales in 
year one. Again we can see the represented effects of natural and fishing 
mortalities. There are two unknowns, F and M, affecting the age structure 
и љенроршћа топ. ethus ,wifesansestamate.can beemade of one, then both in 
combination, the other can be determined. To accomplish this end, the 
ratio of the six year olds and older in year one, to the seven year olds 
and older in year two (term 1 of the basic equation (5)), is compared to 
the ratio of the seven year olds and older to the six year olds and older, 
in the base year where only natural mortality has been affecting the age 
meructure of the population (term 2 of the basic equation (5)). 

Finally, an adjustment must be made to compensate for the variation 
in abundance from year to year. It is obvious that if the same amount of 
erort is expended, the catch PEU Ae oa y car of greater 
abundance will be greater than the catch per unit effort in a year of 
lesser abundance. Therefore, the ratio of the C/E in year one to the 
C/E in year two, provides a measure of the relative abundance in the two 
years (term three of the basic equation (5)). Thus by comparing the 
ratios between successive years with the ratio existing in the base year, 
an estimate of the natural mortality is obtained and, since this is 
assumed to be constant, the total fishing mortality may be obtained. 

ЧО The Effects of Migration 

As a result of the seasonal migration, in the case of baleen 
whales, and the more general spreading in the case of sperm whales, there 


exists a north-south age segregation of the population for a major 
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portion of the year (Ohsumi 1966, Nishiwaki, 1967, Nishiwaki 1966). Older 
males penetrate further into the waters of the North Pacific followed by 
younger males who uSually preceed the population front of the migrating 
main body of whales. There then exists a high probability that an 
examination of the age structure in any ore ten degree square will not be 
representative of the age structure, but will yield a disproportionately 
high number of one particular age group. To eliminate the effects of 
migratory segregation of the population, the age structure and effort are 
considered for ten degree squares in geographic areas. These areas are 
composed of twenty degrees of longitude extending from the equator to the 
pole, Figure 2, thus considering the age structure at all latitudes. 
Then, for example, for area III between 140°W - 160°W: 


= + Е: Х 
m 2773) 28 + + P28 \ 


and 


с ре и 2 
111 5 “m27 * wes * + Сов 


Since the yearly totals of X, C, C/E, etc. are used in the estimation 
of q and these yearly totals are the sum of the individual Squares. au. 
may be applied to the individual squares to yield monthly estimates of 


population. This approach, however, requires extensive hand calculations. 


B. AGE DETERMINATION 

Application of the "д" method requires that the age groups of the 
catch be analyzed with respect to the length-frequency distribution. The 
relationship of length to age was determined indirectly from frequency 
diagrams, Figures 8-19, relating body length to the count of ear-plug 


laminations. 
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Figure 8. Male Fin Whale Body Length vs Count of Ear-Plug Laminations 
(Compiled by Biologists at the Japanese Whaling Institute) 
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Нїриге 9. 
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Male Fin Whale Body Length vs Count of Ear-Plug Laminations 
(Compiled by Biologists at the Japanese Whaling Institute) 
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Figure 10. Male Fin Whale Body Length vs Count of Ear-Plug Laminations 
(Compiled by Biologists at the Japanese Whaling Institute) 
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Figure 11. Male Fin Whale Body Length vs Count of Ear-Plug Laminations 
(Compiled by Biologists at the Japanese Whaling Institute) 
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Figure 12. Female Fin Whale Body Length vs Count of Ear-Plug Laminations 
(Compiled by Biologists at the Japanese Whaling Institute) 
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Erpure 15. 


Female Fin Whale Body Length vs Count of Ear-Plug Laminations 
(Compiled by Biologists at the Japanese Whaling Institute) 


Numbers Indicate Totals of Observations 
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Figure 14. Female Fin Whale Body Length vs Count of Ear-Plug Laminations 
(Compiled by Biologists at the Japanese Whaling Institute) 
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Eppure 15. 


Female Fin Whale Body Length vs Count of Ear-Plug Laminations 
(Compiled by Biologists at the Japanese Whaling Institute) 


Numbers Indicate Totals of Observations 
Area V 
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Figure 16. Male Sei Whales Body Length vs Count of Ear-Plug Laminations 
(Compiled by Biologists at the Japanese Whaling Institute) 
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Figure 17. Male Sei Whale Body Length vs Count of Ear-Plug Laminations 
(Compiled by Biologists at the Japanese Whaling Institute) 
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Figure 18. Female Sei Whale Body Length vs Count of Ear-Plug Laminations 
(Compiled by Biologists at the Japanese Whaling Institute) 


Numbers Indicate Totals of Observations 
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Figure 19. Female Sei Whale Body Length vs Count of Ear-Plug Laminations 
(Compiled by Biologists at the Japanese Whaling Institute) 
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A best fit curve was determined for each length lamination frequency 
diagram by the least squares method using the LSQPL 2 subroutine in the 
W.R. Church Computer Center library. This subroutine has a maximum 
Capacity of one hundred points. For length-frequency diagrams having 
more than this number of observations, the mean value of body length was 
calculated for each number of laminations and the least squares curve 
fitted to these points. 

iiesindividual curves obtained Were then plotted, Figures 20-23, 
and a mean curve determined by applying the known lengths and ages of 
physical and sexual maturity to these curves. The final curve is 
approximately exponential in form and agrees reasonably well with von 


Bertalanffy's theoretical growth equation (Mackintosh 1965, Laws 1959): 


1, = L т кр 
where: 
1, = length at any age 
L = mean length at physical maturity 
k = a constant determined by the physiological characteristics 


of the animal 
A computer program deSignated AGEPROP (see Appendix B) was written 
to calculate the proportion of each age group in the catch, based on 


the age-length curves. 


fee EXPECTED CATCH METHOD 

The basic concept of this method is to obtain an estimate which mini- 
mizes the sum of the squares of the differences between actual and 
expected catches. The unknown parameters are the initial population 


size, N, and the ratio, q, between the instantaneous fishing mortality, 
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F, and the fishing effort, X. Natural mortality is assumed constant, 
and recruitment is obtained by examination of the age structure in 
successive years. The method constructs the population in successive 
years by taking the initial population, subtracting the actual catch, 
subtracting mortality, and adding recruitment. This process is 
repeated in successive years using each estimate obtained as the initial 
populatzonsfor the following year. The expected catches are then 
obtained for each year by dividing the 21 average population 
during the season by the estimated fishing mortality rate. The mortality 
rate is obtained by multiplying q times the known effort. The method 
assumes that catching, natural mortality and recruitment proceed 
successively thereby simplifying the calculations. The essence of the 
method is paraphrased here; for further detail the reader is referred 
to the original paper (Allen, 1966). 

Mathematically, letting C. and X. be catch апа effort respectively, 
pen: 

N, - initial population 

(N,-c) eH - survival to the next year 

(Ni - 21) qX; = expected catch 


Then the following year (№): 


-M 
(N, - суе 


(1 - W) 


where W is the number of recruits in the next year. 
Mis the natural mortality and 


Es 15 the catch in year 1. 
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Survival to the beginning of the next season 


(N -C1)e l m 
Ë (1-W.) Уа | 
2 
and expected catch 
(№1-С1)е  £, 
(1-2) 


Continuing in this manner we сап 


the population equals: 


show that at the beginning of year t 





1 t-1 
С. x (1-W. C x (1- w ) 
б-н о азу ӘЗ) Sey P 
ИЕ Zo Ее 
t t өті E 1)m o7 (t-2)m 
X (1-W.) 
Е 1 
П 
= N -f(C 1 
t шиг 
- (t-1)m 
where: A = # 
t + 
2 (1-9, ;) 
1552) 
ап4 
t-1 C: 
£(C) - С + X 
t-1 1 1-2 A; 


In year t the difference between 


by: 


С N. ot 
NU. - БСХА 


The sum of the squares of the differences, 


year T equals: 


the actual and expected catch is given 


qX+ 


x ах, № 


d, for the period year 1 to 


48 





ах, 2 
) - ша = d 
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ChI +4 
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Expanding and regrouping terms yields: 
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where all terms are summed from t = 1 to t = T. 
The normal equations for which the sum of squares is a minimum are 
obtained by differentiating with respect to N, and q and setting equal 


ШО лего. Thus: 


2 
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Solving for N. and q gives: 
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_ RP - 2SL 


q 
4SU - R? 


Once М, has been estimated, the populations for subsequent years can be 


estimated by: 
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ГІ КЕЗІНЕ 


Figures 24-32 are average monthly estimates of the number of whales 
which will be encountered as false sonar targets in 1000 nautical miles 
of steaming with a 1000 yard sonar range, assuming ideal sonar conditions. 
The values for expected false targets with other sonar ranges may be 
obtained by multiplying these estimates by the sonar range in kiloyards. 

It can be seen that the false target threat exists in three major areas: 
a northern area extending across the North Pacific encompassing the Aleutian 
Islands and the Bering Sea; a central area adjacent to the coast of Japan; 
and a southern area in the Central North Pacific from 30°N-40°N which 
broadens to 50°N in the area east of 160° W. 

. The concentrations of whales as false targets in the northern region 
gradually increases and spreads from west to east beginning in April with 
the northern migration of fin and sei whales along the coast of Japan into 
the Aleutian waters and Bering Sea. The maximum of sixty-three false 
targets per 1000 nmi is reached in the southern portion of this northern 
area in August. At this time, the feeding fin and sei whales are accompa- 
nied by the maximum of male sperm whales which have migrated from the 
equatorial regions. In subsequent months, the number of whales in the 
northern region decreases and shifts westward as the whales migrate south- 
ward to their winter grounds in warmer waters. This area is virtually 
empty of whales from December to late March. 

In the central region adjacent to Japan, the minimum concentration 
also occurs during the period January through March. The population 


bégins to increase in early April due to the northward migration of fin 
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and sei whales and later with the more general spreading of male sperm 
whales. 

The central region exhibits a bimodal distribution; the first peakin June 
is due to the northern migration of fin and sei whales and the second 
maximum is in October and corresponds to the peak southern migration. 
This second maximum of thirty whales per 1000 nmi is greater than the 
first of 9 per 1000 nmi. This is because the southward migration of both 
fin and sei whales closely coincides while the northern migration is 
distributed over a longer period of time (Kellog 1929). 

The southern area shows its initial increase in May in the eastern 
region due to the arrival of sperm whales along with the appearance of 
Some tin whales off the west coast of the United States. The population 
of this region reaches a maximum in July, coincident with the passage of 
fin whales into ie northern region. Concurrently, young male, female 
and calve sperm whales arrive in the Central North Pacific and spread 
along 30°N. The sperm whale maximum occurs in Spetember and then de- 
creases until, in November, this area is also virtually empty of whales. 

The total estimated populations of fin, sei and sperm whales in the 
North Pacific north of 30°N, for the years 1966-1970, are listed in 
Table II. Results are provided for both the "q" and Expected Catch 
methods of population dynamics. Areas III through VI contain estimates 
based on combined Japanese and Russian data, while Area II contains 
Russian data only. The blank squares represent areas for which the 
particular method produced anomalous negative results and do not appear 
in the table. Direct comparison, from the table, of the results gener- 
ated by the two methods is difficult since the Expected Catch method 


yields the population at the start of the season and the ''q' method 


yields the population at mid-season. 
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TABLE III 


CALCULATED POPULATIONS OF FIN, SEI AND SPERM WHALES IN THE 
NORTH PACIFIC BY THE "а" AND EXPECTED CATCH METHODS 


1966-1970 
AREA II - V 


Calculated Populations 


Ba Expected Expected Expected 
Па нов Catch Mond Catch Method Catch 





AREA II 
8423 1124 59515 
8813 1120 40399 
7188 1079 40461 
4113 710 38939 
5265 613 5/191 
Average population N; 
6761 ЕЭ 59521 
AREA III 
12402 1299 
4890 449 
5038 360 
1948 264 17500 
4026 211 
Average population N; 
5661 517 17500 
AREA IV 
5589 2719 
5526 1966 
9856 1677 22100 
5662 1525 27300 
1991 1181 
Average population N; 
4521 9068 24700 
AREA V 
21257 4853 20177 
5185 4554 20724 
4302 4034 21416 
4075 5614 21970 
2252 5460 22079 
Average population N; 
3714 4103 21-12 
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TABLE IV 


CALCULATED POPULATIONS OF FIN, SEI AND SPERM WHALES IN THE 
NORTH PACIFIC BY THE 'q'' AND EXPECTED CATCH METHODS 


1966-1970 


AREA VI 


Calculated Populations 


ша)! | Expected Expected Expected 
Method Catch Е Catch eed Catch 


AREA VI 


2221 1530 
6340 1506 
5856 LOSS 
10382 863 
4448 835 


Average population N; 
6049 1151 
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IV. DISCUSSION 


A. DATA 

The data used in this thesis were extracted from the basic raw data 
reports prepared by Japanese and Russian whaling institutes for inclusion 
те International Whaling Statistics. In their basic form, these 
Reports contain, for each ten degree square, M listings of catch 
by species, sex and catcher effort expended. An annual length-frequency 
distribution by statistical area is also included. By using these 
reports, information may be derived for monthly population estimates. 
Whaling data presented in the International Whaling Statistics are or- 
ganized by geographic area and seasonal totals and can only be used for 
estimates of total annual population of the North Pacific. 

A serious limitation exists in using whaling data for monthly popula- 
tion estimates in that not all areas are fished in every month. The 
whaling season is limited by regulation of the International Whaling 
Commission to six months annually for baleen, and eight months for sperm 
whales. This limitation is compensated for in part in that the whaling 
fleet follows the migration of the whale herds and seeks out areas of 
Significant concentration. 

A second problem which exists lies in the minimum size limitations of 
eiescommission which are 57 feet for fin (Balaenoptera physalus), sei 
(Balaenoptera borealis), and sperm (Physeter catadon) whales. Thus the 
length frequency distribution may not truly reflect the age composition. 
of the catch since certain ages, determined by length, are excluded from 
Euescestchable population. This may bias the estimates of recruitment used 


in the dynamical methods. 
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Ihird, due to the End nature of the fishery, data relating to 
the abundance of whales outside the normal season is virtually non-existent. 
This results in a period of a year for which no population estimates may be 
calculated. In this study, this period comprises the months of January, 
February and March. However, examination of the migratory cycles of the 
major whale types (Fig 33, Kellog 1929) indicates that very few whales 
are found in these northern waters in the winter months. The few that 
persist consist of old sperm and fin bulls which can withstand large 
temperature extremes. As these whales are usually loners, they should 
Present no significant false target threat. 

цаша the data represents only the catch statistics for the 
Mier ciatly exploited species, fin, sei, bryde's, sperm and, for earlier 
years, humpback whales. The ramifications of this problem are discussed 


further below. 


B. AGE DETERMINATION 

Age determination is important for development of accurate methods in 
the fisheries biology of whales, and has been the subject of intensive 
investigation. Of all the age characteristics studied, the accumulated 
laminae in the ear plug of the external auditory meatus is believed to be 
most valuable criteria for age determination of the catch (Ichihara 1966). 
The conical shaped plug consists of alternating layers of dark keratinized 
cells and bright degenerated fatty cells which form at a rate of one to 
eyo per year. It is believed that the formation of alternating layers is 
in response to the biannual hormonal cycle and periodic changes of food 
supply. It is assumed that vitamin A, plentiful food supply, and growth 
hormones result in the formation of a bright layer. Vitamin A deficiency, 


fasting and estrogen result in formation of a dark keratinized layer. 
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The bright layer in the sexually mature whale is formed during the period 
of feeding migration and the dark layer during a breeding migration. 

The annual rate of formation of ear-plug laminations has been a source 
of disagreement among biologists. Several of the possibilities offered 
are: 

1) one lamina formed per year 

2) two lamina formed per year 

3) two lamina formed per year until sexual maturity is reached and 
one lamina per year thereafter 

4) 1.5 lamina formed per year. 

For the calculations carried out in this paper, age structure was based 
upon an assumed rate of lamina formation of 1.5 per year for the sexually 
immature whales. This figure was selected as the most probable estimate 
in view of work done 20202 (1966) who verified it by three methods 
and Oshumi (1964) who reached the same conclusion by investigation of 
natural mortality rates. As the same standard is used for age determina- 
tion of catches in successive years for comparison of year groups, any 
error introduced is constant in the determined age structure, and is 
insignificant in the final calculations. Allen (1967) calculated the 

fin whale population of the North Pacific using three conversions from 
ear-plug laminations to years (one per year, one per year to five years 
and then two per year and two per year). His results showed an average 


difference of only two percent between the calculated populations. 


C. AREAS OF DISTRIBUTION 
A significant point in the consideration of whale distributions is 


ШЕ identification of the specific areas of distribution. These areas 


> 
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must be arrived at by examination of the catch statistics in conjunction 
with the migratory cycles of the animals of interest. Whales migrate in 
response to two factors: 

1) location of areas containing specific foods in sufficient 
concentration 

2) the requirement for a suitable environment for the production and 
rearing of young. 
штап бе seen that the composite areas of distribution, Figure 6, closely 
agree with the zooplankton distribution, Figure 34. This is expected 
since it is known that the abundance of food controls the limit of migra- 
tion (Nemoto 1959). 

The oceanic structure is intimately associated with the production of 
zooplankton and, therefore, determines the location of whaling grounds. 
Vertical and horizontal temperature gradients and vertical stability of 
sea water masses play an important role in the replenishment of micro- 
nutrients, mainly phosporic and nitrogen salts. These salts are the 
limiting factors in the production of phytoplankton, the food supply of 
zooplankton. Therefore, the highest concentration of whales, and thus 
the most troublesome false target areas, are those where cold and — 
water masses meet and the three basic conditions of light, nutrients and 
oxygen are present (Ruud 1938, Nasu 1963). It should be noted that whale 
schools do not concentrate in the rapid, easterly rapidly flowing zones 
of warm water. Therefore, few whales are found in the west wind drift. 

The accuracy of the areas of whale distribution (Figure 6) is 
substantiated by the fact that in the 1969 season, only 183 whales out of 
a total catch of approximately ten thousand were taken outside these areas 
in spite of the efforts of scouting boats to locate new whaling grounds 


(International Whaling Commission 1969). 
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D. SONAR PROBABILITY FACTOR 

Ideal sonar conditions were assumed in this study; therefore, a value 
of 1 was used for the sonar probability factor in the application of the 
modified "Transect Method.'' Since the effective sonar range is a function 
of the thermal structure, density and salinity characteristics of the 
oceanic region, the sonar probability factor will vary seasonally with 
the meteorological and oceanographic conditions. 

Bow unitS Operating in the North Pacific provide the best immediate 
source of data to evaluate the sonar probability factor and accurately 
determine the values associated with specific oceanic regions. A feed- 
back system from operational units would expedite the resolution of this 
problem. 

It is known that the migration and feeding grounds of fin and sei 
Whales are a function of the sea surface temperature and the locations 
of water mass fronts (Uda 1954, Omura 1955, Uda 1956, Nasu 1957, Uda 
me? Uda 1958, Nemoto 1959, Nasu 1960, Uda 1962, Nasu 1966). It may, 
mieretore, be possible to determine the sonar probability factor through 
Fleet Numerical Weather Central's thermal structure analysis program. 

Use of this program may provide prediction of the location of the main 
body front of the migrating herds of whales, and determination of favorable 
feeding grounds. 

The prototype atlas produced by Leapley and Leavenson (1969) bases 
prediction of false targets on sighting reports by aircraft and shipping. 
This approach is limited by several factors. First, whale sighting, 
especially from aircraft, is seriously hampered by adverse meteorlogical 
conditions which reduce visibility. Second, the characteristic diving 
behavior of whales may result in their being unobserved from a rapidly 


moving platform, i.e. a survey aircraft. These factors could result in 
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a lesser number of whales being observed in an area than are actually 
present. On the other hand, when whales are sounding and surfacing the 
same whale may be counted more than once. This would result in an 
erroneously high count for the area under observation. Finally, this 
approach is limited in that no provision is incorporated for variations 
EnESonar conditions. It is felt that the approach taken in this thesis, 
the analysis of whaling data, can provide more accurate estimates of the 
false target threat if it is combined with oceanographic research to 
determine the value of the sonar probability factor. Both can be sent 
puro the Fleet as a part of the Fleet Numerical Weather Central 


forecast. 


P INCLUSION OF OTHER WHALE SPECIES 

This thesis considers only the species of whales hunted and reported 
on by the commercial fisheries, i.e. fin, sei, brydes (included under 
sei) and sperm whales. Other whales of importance to the false target 
problem in the North Pacific include: the Northern Right (Eubalaena 
sieboldii), Gray (Eschrichtius gibbosus), Pigmy Sperm (Kogia breviceps) 
and Humpback (Megaptera novae-angliae). Further study and analysis of 
the Weer ibution and migratory patterns of these whales would provide 
further refinement of the analysis method presented in this paper. 

For example, a census is being conducted this year of the gray whale. 
At Point Yankee, California, the southern migration reached a peak with 
the passage of two hundred whales per day (Dr. Dale W. Rice, private 
communication). The daily count of the passage of the whales past a 
Specific point such as Point Yankee, would yield the distribution of the 
whale population in migration. This distribution may be combined with 


the known migratory swimming speed of four knots and the known migratory 
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Paths along the coast of North America into the Bering Sea (Pike 1962, 
Rice 1971). Thus an estimate of the arrival times of the main body in 


specific areas can be derived. 


Е. WHALE SIGHTING PROGRAM 

In the Fall of 1971, Dr. Paul Sund of the National Marine Fisheries, 
ШЕ tiated a whale sighting program for the Pacific Ocean (Dr. C. Fiscus, 
Marine Mammal Laboratory, Seattle, Washington, personal communication). 
Under this program, officers of United States merchant ships are being 
moment to recognize species of whales 228. through the water and 
report whales sighted during transits of the Pacific Ocean. These sighting 
reports will be card punched for computer processing and will provide a 
valuable source of data on whale distribution patterns when a sufficient 


quantity of data has been accumulated to give meaningful results. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The analysis of commercial whaling data for the North Pacific provides 
a useful method for estimating the false target threat in this area. 

Studies should be undertaken to determine the values of the sonar 
probability factor for those areas identified with significant whale 
density. This may be accomplished by means of feedback reports from 
operational units combined with Fleet Numerical Weather Central's thermal 
structure forecast. These studies should also incorporate the results of 
the new whale sighting program initiated by Dr. Paul Sund of the National 


Marine Fisheries, in 1971. 
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